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ABSTRACT 

We report the results of the hrst X-ray observation of the luminous and helium-rich O-type sub- 
dwarf BD+37°442 carried out with the XMM-Newton satellite in August 2011. X-ray emission is 
detected with a flux of about 3xl0 -14 erg cm~ 2 s _1 (0.2-1 keV) and a very soft spectrum, well 
fit by the sum of a blackbody with temperature /cTbb = 45^ 9 eV, and a power law with a poorly 
constrained photon index. Significant pulsations with a period of 19.2 s are detected, indicating that 
the X-ray emission originates in a white dwarf or neutron star companion, most likely powered by 
accretion from the wind of BD +37° 442. 

Subject headings: Stars: individual: BD +37° 442 — subdwarfs — pulsars: general — X-rays: stars — 
X-rays: binaries 



1. INTRODUCTION 

Hot subdwarfs are evolved, low-mass stars that have 
lost most of their hydrogen envelop e and are now in the 
stage of helium-core burning (see iHeberl 120091 for a 
review). From the spectroscopic point of view we dis- 
tinguish the cooler B-type subdwarf (sdB) stars, with 
effective temperature T e ff < 40,000 K, an d the hotter 
O-ty pe (sdO) stars, with T cS > 40,000 K (jHirsch et al.l 
2008) . Most sdB stars are helium poor and display no or 
only weak helium lines, while most sdO stars are helium 
rich and show He I and He n lines. While the sdBs form 
a homogenous group, the sdOs show a wide spread in 
temperature (T eff = 40,000-100, 000 K), gravity (}or(q) 
= 4-6.5) and heli um abundance (|Heber fc Jeffervl [19921 : 
IHeber et al.l 120061 ). Historically, sdO stars were divided 
into 'luminous' and 'compact' sdOs, depending on their 
value of log(g). 

One possible mechanism responsible for the loss of 
the massive hydrogen envelopes of hot subdwarfs is 
mass transfer in a binary. Indeed, there is substan- 
tial evidence th at many hot subdwarfs are in close 
binary systems dMaxted et al.l 120011 : iNapiwotzk i et al.l 
120041 : iCopperwheat et al.l 1201 lh . In particular, a large 
fraction of binaries is found among sdB stars. Models 
of binary evolution ([Han et al.l 120021 ) predict that many 
of the subdwarf companions should be w hite dwarfs 
(WD ). Indeed a few sdB + WD are know n dKpen et all 
[llMlMaxted et al.ll2000t lEdelmann et alll2005l )." and re- 
cent optical studies of single-lin ed spectroscopic b inaries 
yielded several new candidates ([Geier et all 12010). The 
presence of a compact companion (a WD, neutron star 
(NS), or black hole) could be revealed by the detection 
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of X-ray emission powered by accretion, if the subd- 
warf mass-donor can provide a sufficiently high accre- 
tion rate. Short Swift/XRT observations of candidate 
sdB+W D/NS binaries, selected from the iGeier et all 
(2010) sample, gave X-ray luminosity upper limits of 
- 10 30 - 10 31 erg s" 1 (jMereghetti et all 1201 lah . These 
limits confirm that sdB stars have rather weak stellar 
winds, unable to provide enough accretion rate. 

Although the fraction of binaries among sdO stars is 
not as high as for the sdBs, the prospects to find X- 
ray emitting companions are more promising for these 
stars. In fact, at least a few luminous sdO stars show 
evidence for stellar winds with mass loss ra te M ~ 
10- 7 -10- 10 M o yr" 1 (jJefferv fc Hamannll2010l ) . Among 
these, the bright star HD 49798 is known as a soft X- 
ray source, with a strong periodic modula tion at 13.2 s, 
since more than 15 years ([Israel et al.lll997[ ). More recent 
X-ray observations with XMM-Newton showed that the 
compact companion of HP 49798 is a massive (1.28 ± 
0.05 M Q ) white dwarf (jMereghetti et alll2009D . To our 
knowledge, no other X-ray detections of sdO stars have 
been reported up to now. 

Prompted by our findings on HD 49798, we carried out 
an X-ray observation of anothe r brig ht sdO star. Our 
target is BD +37° 442 (|Rebeirotj[l966( ). a luminous sdO 
which is often referred to as an 'extreme helium star' be- 
cause of its complete lack of hydrogen (while HD 49798 
has about 20 % (by mass) of hydro gen): it shows ev- 
idenc e for significant mass loss rate (jJefferv fe Hamannl 
I201CO and stellar parameters similar to those of HD 49798 
(see Tab. dj, but with the notable difference that, up to 
now, no evidence of binary nature was reported. Here 
we present the results of the first X-ray observations of 
BD +37° 442, which indicate the presence of a compact 
companion, either a WD or a NS. 
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TABLE 1 

Main parameters of the sdO stars BD +37° 442 and HD 49798. 



Parameter 


Symbol 


BD +37° 442 


HD 49798 






Value 


Reference 


Value 


Reference 



Mass (M ) M 0.9 

Radius (Jig) R 1.6 

Surface gravity log g 4.0 ± 0.25 

Luminosity (L ) L 25,000 

Effective Temperature (K) T cff 48,000 K 

Magnitudes U, B, V 8.57, 9.73, 10.01 

Distance (kpc) d 2+° '% 

Terminal wind velocity (km s _1 ) Voo 2,000 

Mass-loss rate (M Q yr" 1 ) M 10 



-8.5 



1.50 ± 0.05 
1.45 ± 0.25 
4.35 
14,000 
46,500 
6.76, 8.02, 8.29 
0.65 ± 0.1 
1,350 

10 -8.5 
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2. OBSERVATIONS AND DATA ANALYSIS 

BD +37° 442 was observed with XMM-Newton on 2011 
August 21, starting at 23:35:19 UT (MJ D = 55794.983) 
The three EPI C cameras, i.e. one p n intruder et al1l2001[) 
and two MOS ([Turner et al.ll200ll ). were operated in full 
frame mode, with time resolution of 73 ms for the pn and 
of 2.6 s for the two MOS cameras; the effective source 
exposure time was, respectively, of ~ 33 ks and ~ 28 ks. 
For all cameras the medium thickness filter was used. 

We used version 11.0 of the XMM-Newton Science 
Analysis System (SAS) to process the event files. The ob- 
servation was affected by high instrumental background 
only for a short time interval of ~ 1.2 ks, which we re- 
moved for the spectral analysis. We selected only the 
events with pattern in the range 0-4 (i.e. mono- and bi- 
pixel events) for the pn camera and 0-12 (i.e. from 1 to 
4 pixel events) for the two MOS. 

The images obtained with the three cameras clearly 
show that BD +37° 442 is emitting in the X-ray range: in 
fact a source is significantly detected at the coordinates 
R.A. = 01 h 58 m 33.4 s , Dec. = +38° 34' 22.0", which dif- 
fer by only 1.8" from the position of BD +37° 442. This 
difference is consistent with the ~ 2" r.m.s. astromet- 
ric accuracy of XMM-Newton Q. The source net count 
rate in the 0.15-2 keV range is (1.62 ± 0.08) xlO -2 cts 
s _1 and (1.9 ± 0.2) xl0~ 3 cts s _1 for the pn and each of 
the two MOS, respectively. BD +37° 442 is not detected 
above ~ 2 keV, suggesting a soft spectrum, as confirmed 
by the spectral analysis presented below. 

For the timing analysis, we used the data of the whole 
observation and the three EPIC cameras, extracting the 
events from a circular region with radius 15" centered 
at the source position. We converted the arrival times 
to the solar system barycenter, then we combined the 
three datasets in a single event list. The background- 
subtracted light curve of BD +37° 442 does not show 
variability on time scales from hundreds of seconds to 
the observation length. In order to look for periodicities, 
we examined 65,232 statistically independent periods in 
the r ange 1-10,000 s app lying the phase-distribution 
test (|Buccheri et al.l fT983r> to the counts in the 0.15-2 
keV energy range (685 counts). The resulting distri- 
bution of Z 2 statistics versus trial period is shown in 




0018.ps. gz 



http: // xmm2 . esac.esa. int /docs/ documents /C AL-TN- 



FlG. 1. — Distribution of the Z 2 statistics versus trial period for 
the BD +37° 442 events in the 0.15-2 keV energy range. 

Fig [TJ A peak value of Z 2 = 35.15 is found at the 
period P = 19.16 s. The corresponding probability of 
chance occurrence, after taking into account the number 
of trials, is 1.5 x 10~ 3 . As shown by the folded and 
background-subtracted light curves of Fig. [2j the peri- 
odicity is visible also in the separate data sets, i.e. pn 
and MOS1+MOS2, which give consistent phase distri- 
butions and pulsed fractions. The refined period value, 
obtained through pulse phase fitting, is P = 19.156 ± 
0.001 s. The pulse profile shows a single broad peak 
with a 0.15-2 keV pulsed fraction of 31 ± 4 % (ampli- 
tude of a sinusoid fitted to the folded light curve divided 
by the average flux). No significant differences are seen 
in the light curve shape and pulsed fraction in differ- 
ent energy ranges. We also performed a phase-resolved 
spectral analysis, as described below, for pairs of spec- 
tra corresponding to the pulse maximum/minimum or 
rise/decrease phase intervals, without finding significant 
differences in the spectral parameters. 

Due to the limited signal-to-noise ratio of the MOS 
data, for the spectral analysis we considered only the pn 
events. We generated the applicable response matrices 
and ancillary files using the SAS tasks rmfgen and arf- 
gen. To ensure the applicability of the % 2 statistics, the 
spectrum was rebinned with a minimum of 30 counts per 
bin and fitted using XSPEC (V 12.7.0). We only used 
the energy range 0.2-1 keV since above 1 keV the back- 
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luminosity is obtained for the lowest temperature). The 
0.2-10 keV luminosity of the power-law component is in 
the range (1.8-3.2) x 10 31 erg s _1 . 



Phase 

Fig. 2. — Background-subtracted light curves of BD+37° 442 in 
the energy range 0.15-2 keV, folded at the best-fit period P = 
19.156 s. From top to bottom: pn data, sum of MOS1 and MOS2 
data, total data (pn+MOSl+MOS2). 



ground dominates and the source flux is negligible. In 
the following, all the spectral uncertainties and upper 
limits are given at the 90 % confidence level for one in- 
teresti ng parameter, and we assume a source distance of 
2 kpc (IHusfeldlll987h . 

The source spectrum is very soft: a fit with an ab- 
sorbed power law gives an unphysically large photon 
index (r ~ 5), while a blackbody gives a temperature 
feTsB = 93 eV. However, these single-component mod- 
els give unacceptable fits (x 2 > 2). A model consist- 
ing of an absorbed power law plus blackbody provides 
a significant improvement (xt ~ 1-24 for 9 degrees of 
freedom, Fig|3]), although the power-law photon index 
r is poorly constrained in the range between -1 and 4. 
Keeping T fixed at its best fit value (2.25), the values 
of the other spectral parameters are: column density 
Nh = (lAto'l) x 10 21 cm" 2 , blackbody temperature 



on+162 
J3 -28 



kTsB = 45_ 9 eV, and emitting radius Rbb 
km. The absorbed flux in the energy range 0.2-1 keV is 
/abs,x = (2.6 ±0.3) x 10- 14 erg cm" 2 s" 1 , of which 28 % 
can be ascribed to the power-law component. The bolo- 
metric flux of the blackbody component, corrected for 
the absorption, is 1.7 x 10 -12 erg cm -2 s _1 , correspond- 
ing to a luminosity of 7.8x10 erg s" 1 (for a source 
distance d — 2 kpc). This value is subject to a large 
uncertainty because the value of /cTbb is only poorly 
constrained and depends on the slope of the power-law 
component. Therefore, in order to assess in a conser- 
vative way the uncertainty on the source luminosity, we 
repeated the spectral fits fixing the photon index at the 
values r = 1 and F" = 3. In this way, we found accept- 
able values of /cTbb in the range 30-68 eV (90 % c.L), 
corresponding to bolometric blackbody luminosities be- 
tween 6xl0 31 and 1.5xl0 35 erg s _1 (where the highest 




Fig. 3. — Top panel: pn spectrum of BD +37° 442 with the best- 
fit power-law (dashed line) plus blackbody (dotted line) model. 
Bottom panel: residuals (in units of a) between data and model. 



3. DISCUSSION 

Our XMM-Newton observation of BD+37°442 has 
provided the first detection of X-ray emission from this 
O-type subdwarf. Shocks in the radiation-driven stellar 
winds of main sequence, giant, and supergiant early type 
stars give rise to soft X-rays, with a typical ratio of X-ray 
to optic al flux of th e order of log(f!}%_? kcV //boi) = -6.45 
± 0.51 (Naze 2009). It is not known if this relation holds 
also for low luminosity stars. The only other sdO star 
possibly detected in X-rays is HD 49798, which is in a bi- 
nary system with an X-ray emitting white dwarf. In this 
system, the X-ray flux of fx — 4.3 x 10~ 14 erg cm -2 s" 1 , 
visible when the accreting WD companion is in eclipse, 
corresponds to log(/^ ;|^ keV //boi) = -7.4 and could be 
due to emission from the sdO (|Mereghetti et al.ll2011b| ). 
For BD +37° 442 we measure log(/ u ^ b 1 s keV //boi) = - 
6.63, consistent with the above values, but the signif- 
icant periodic modulation at 19.156 s clearly indicates 
that most of, or possibly all, the detected flux originates 
in a compact companion star, rather than in BD +37° 442 
itself. The periodic modulation can be explained equally 
well as the spin period of a NS or of a WD. 

Ultraviolet spectra of BD +37° 442, acquired with the 
IUE satellite, show Nv and Civ resonance lines with P 
Cygni -like profiles, ind icating the presence of a stellar 
wind (|Rossi et al.l 11984). A more recent analysis of UV 
an d optical spectra o f BP + 37° 442 has been reported 
by Uefferv fc Hamannl (|2010l ). who derived a mass loss 
rate of 3xl0 -9 M yr _1 and a wind terminal velocity 
of Voo = 2,000 km s _1 . It is possible that part of the 
sdO stellar wind is captured by its compact companion, 
giving rise to accretion powered X-ray emission. The 
density and velocity of the wind material at the posi- 
tion of the compact object can be computed assuming 
a canonical wind velocity law with radial dependence 
v(R) = Voo(l — 1/R)P, where R is the radial distance 
in units of stellar radii and the index /3 is typically in the 
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range 0.6-1. Assuming for simplicity Bondi-Hoyle accre- 
tion, and relating the orbital separation to the period of 
a circular orbit P or b through Kepler's law, we can derive 
the expected accretion luminosity as a function of P or b 
shown in Fig. 01 The curves refer to a NS of 1.4 M Q and 
10 km radius, and to a WD of 0.6 M and 10,000 km 
radius; i n both cases w e used for BD +37° 442 a mass of 
0.9 M Q (|HusfeldlH987l ). We have also indicated in the 
figure the allowed range of luminosity of the blackbody 
component (for simplicity we do not consider the addi- 
tional 30 % contribution to the total luminosity coming 
from the power-law component). 



10' 



spectroscopic (iFay et al j|1973t iKa ufmann fc Thelllll980fc 
Dworetskv et all I1982D or photometric data ( Landold 
19681 [1973T ). and infrared observat ions did not show any 
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Fig. 4. — Estimated wind-accretion X-ray luminosity of 
BD +37° 442 as a function of the binary period, in the case of 
a WD (dashed lines) and of a NS (solid lines) The couples of lines 
refer to two values of the wind velocity index (3 = 0.6 and /3 = 1. 
The horizontal dotted lines represent the minimum and maximum 
of the estimated luminosity of the blackbody component. 

Although the observed luminosity is only poorly con- 
strained due to the large uncertainties in the spectral 
parameters, it is clearly consistent with that expected 
from a NS orbiting BD +37° 442 with a period between 
about one and several days. The blackbody emitting ra- 
dius derived from the best fit, iZeB = 39l2,g 2 km (for d = 
2 kpc), is only marginally consistent with a neutron star. 
However, this parameter is strongly correlated with the 
poorly constrained slope of the power law spectral com- 
ponent and an acceptable fit can be obtained, e.g., with 
T = 1.5, kT-QB ~ 58 eV, and i?BB ~ 10 km. Alternatively, 
the accreting companion could be a WD, if the sdO ex- 
tends to (or close to) the Roche-lobe, thus yielding an 
accretion rate larger than what we computed assuming 
wind accretion (or in the unlikely case that the adopted 
distance is largely overestimated). The observed spec- 
trum is much softer than the typical spectrum of neutron 
stars in classical X-ray binaries, which however have very 
different companion stars and accrete at higher rates. It 
should also be considered that we can observe the low 
temperature emission from BD +37° 442 (and HD 49798, 
see below) only because these sources are at high Galac- 
tic latitude, where interstellar X-ray absorption is very 
small. 

Up to now BD +37° 442 was believed to be a sin- 
gle star: no evidence for a companion was seen in 



excess emission (jTheill et al.lll995| ). On the other hand, 
some of the very few published radial velocity measure- 
ments are inconsistent with a single value, hinting to a 
possible binary motion. A radial velocity V Y = -156.4 
± 1.1 km s -1 is give n in the original discov ery paper 
(|Rebeirotlll966h . while IDrilling fc Heberl (fl987h found V T 
= -94 ± 1 km s _1 . New accurate radial measurements 
are needed to determine the parameters of this binary. 

The X-ray properties of BD +37° 442 are very simi- 
lar to those of the only other X-ray source associated 
with an sdO star: HD 49798. In this binary a mas- 
sive white dwarf accretes matter from the stellar wind of 
its subdwarf companion (|Mereghetti et all [20091 ) . Also 
in the case of HD 49798 the X-ray emission is pulsed, 
with a period (13.2 s) similar to that we have discov- 
ered in BD+37°442, and the soft spectrum is well de- 
scribed by the sum of a black body with /cTrr = 39 eV 
and a power law with T = 1.6 (|Mereghetti et al.ll2011bh . 
HD 49798 has an orbital period of 1.55 days, first discov- 
ered through ra dial velocity mea surements of the sdO 
optical emission (jThackeravl 119701 ) . and an X-ray lumi- 
nosity of 10 32 erg s _1 (at d = 650 pc). 

The fact that the only two sdO stars observed with 
sensitive X-ray instruments turned out to be binaries 
with a compact companion is noteworthy. Although it 
is premature to draw conclusions from such small num- 
bers, this might suggest that the fraction of sdO binaries 
is larger than currently believed on the basis of optical 
observations. Note that, contrary to late type compan- 
ions, which can give a detectable contribution in the spec- 
tral and photometric data, WDs and NSs are too faint 
and completely outshined in the optical/UV by the sdO 
emission. Both BD +37° 442 and HD 49798 belong to the 
subclass of luminous, He-rich sdOs, whose origin and evo- 
lutionary link with other classes of sta rs is still unclear 
(|Napiwotzk i 2008; Justham et al.ll2011[ ). In this respect, 
X-ray observations might provide important information 
to complement optical/UV data. Besides the potential 
of discovering other sdO binaries through the detection 
of pulsations, the high sensitivity of current satellites like 
XMM-Newton and Chandra can probe the X-ray emis- 
sion of single sdOs, and test whether their wind emission 
scales as in early type stars of higher luminosity. 
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